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Directed flow of deuterons, tritons, ^He, and ^He is studied in Au+Au 
collisions at a beam momentum of 10.8 A GeV/c. Flow of all particles is 
analyzed as a function of transverse momentum for different centralities of the 
collision. The directed flow signal, vi{pt), is found to increase with particle 
mass. This mass dependence is strongest in the projectile rapidity region. 
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I. INTRODUCTION 



Anisotropies in the azimuthal distribution of particles, also called anisotropic (directed, 
elliptic, etc.) transverse flow play an important role in high energy nuclear collisions 
In particular, flow of composite particles such as light nuclei is of a large interest for the 
understanding of the nuclear collision dynamics At lower beam energies (kinetic en- 

ergies from about 0.2 to 1.15 GeV per nucleon) flow of light fragments has been studied 
intensively both experimentally |^^1^ and theoretically P-|6|JT5[| (and references therein). 
A theoretically predicted increase of {px)/A, the mean transverse momentum per fragment 
nucleon projected onto the reaction plane, with particle mass was later conflrmed experimen- 
tally [[7|-[T^. Initially, the theoretical prediction of the effect was based on the observation 
that the collective motion of the higher mass fragments should be less sensitive to thermal 
distortions. Later, calculations showed the same effect in models based on a picture 
where the production of light nuclei takes place by the coalescence of nucleons close to 
each other in momentum and conflguration space. Experimental data [TT| for particles with 
transverse momentum greater than 0.2 GeV/c per fragment nucleon were found to be con- 
sistent with momentum space coalescence, but a complete understanding of the effect is still 
lacking. Recently, a measurement of light fragment directed flow in Au+Au collisions at 



AGS energies has been reported by the E802 Collaboration |T6[. Their conclusion was that, 
in the target (pseudo)rapidity region, the maximal (as a function of centrality) values of 
{Px)/A do not depend on the particle species implying a common directed flow velocity. 
In the current paper, we go beyond our previous measurements of anisotropic flow at 



AGS energies [^^-[T9| and present results of the analysis of directed flow of deuterons, tritons, 
^He, and ^He, detected in the E877 spectrometer in Au + Au collisions at a beam momentum 
of 10. GeV/c. Directed flow of all particles is analyzed as a function of the transverse 
momentum for different centralities of the collision. Triton, ^He, and ^He flow patterns are 
analyzed in the beam rapidity region, where these particles can be easily identified. Deuteron 
flow is analyzed also as a function of rapidity, in the rapidity region of 2.2< y <3.4. The 



results of a very similar analysis of proton and charged pion flow were reported in [T^ and 
many technical details relevant for the current analysis can be found there. 
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II. E877 APPARATUS AND FLOW ANALYSIS 



The E877 apparatus is shown in Fig. 1. In the E877 setup, charged particles, emitted 
in the forward direction and passing though a colhmator (—134 < Ohorizontai < 16 mrad, 
— 11 < O^erticai < H uirad), are analyzed by a high resolution magnetic spectrometer. 
The spectrometer acceptance covers mostly the forward rapidity region. The momentum 
of each particle is measured using two drift chambers (DC2 and DC3, position resolution 
about 300 yum) whose pattern recognition capability is aided by four multi-wire proportional 
chambers (MWPC). The average momentum resolution is Ap/p ^3% limited by multiple 
scattering. A time-of-flight hodoscope (TOFU) located directly behind the tracking cham- 
bers provides the time-of-flight with an average resolution of 85 ps. Energy loss in TOFU 
and in a Forward Scintillator array located approximately 30 m downstream of the target is 
used to determine the particle charge. 

The particle identification is performed by combining measurements of momentum, ve- 
locity, and charge of the particle. In the present sample of light nuclei we estimate an 
admixture of other particle species to be no more than 15-20%. Given that particles of the 
admixture (mostly other light nuclei, and some protons in the deuteron sample) exhibit a 
similar flow signal (see below), we arrive at relative errors in the final results due to particle 
misidentification of less than or about 5%. 

The determination of the centrality of the collision and of the reaction plane orientation 
are made using the transverse energy flow measured in the target calorimeter (TCal), and 
participant calorimeter (PCal). Both calorimeters have 27r azimuthal coverage and, com- 
bined together, provide nearly complete polar angle coverage: TCal and PCal cover the 
pseudorapidity regions —0.5 < r] < 0.8 and 0.8 < t] < 4.2, respectively [|T^. 

The azimuthal anisotropy of particle production is studied by means of Fourier analysis 
of azimuthal distributions p0| , ^ , p!7| -[l9[| . This yields the rapidity, transverse momentum, 
and centrality dependence of the Fourier coefficients f „ (amplitude of n-th harmonic) in the 
decomposition: 

d^N d^N 1 d^N ^ o 

d-^p ptdptdydcp Iti ptdptdy 

where the azimuthal angle is measured with respect to the (true) reaction plane. 



Similarly to the analysis presented in |T8|,|T9|, the reaction plane angle is determined 



in four non-overlapping pseudorapidity windows. The 'reaction plane resolution', i.e. the 
accuracy with which the reaction plane orientation is determined, is evaluated by studying 
the correlation between flow angles determined in different windows. Finally, the flow signals 
are corrected for the flnite reaction plane resolution. Details of this procedure are described 
in 



18,19 



III. RESULTS. DIRECTED FLOW OF LIGHT NUCLEI 

We present our results in Figures § and ^ for four centrality regions, selected in ac- 
cordance with transverse energy Et measured in PCal and corresponding to the values of 
(ytop{ET)/(Jgeo ^ 23-13%, 13-9%, 9-4%, and <4% (see Fig. 4 in [|1|]). The value of crtop(Er) 
is obtained by an integration of da/dET from a given value of Et to the maximal one ob- 
served, and the geometric cross section is deflned as (Jgeo = '^'^'q 

(^1/3 + ^1/3)2 _ g_^3 

A = 197 and ro=1.2 fm. 

The amplitude of the flrst harmonic in the Fourier decomposition of the azimuthal dis- 
tribution, Vi{pt), for deuterons is presented in Fig. For comparison the results for protons 



from [19| are shown by open symbols. The deuteron directed flow signal {vi{pt)) is system- 
atically larger than that of protons, especially at rapidities close to the projectile rapidity. 
The relative difference in deuteron and proton flow increases with centrality. The centrality 
dependence of deuteron directed flow depends on rapidity. While at rapidities y < 2.4 flow 
almost disappears at high centralities, in the beam rapidity region it decreases very little. 

At rapidities close to beam rapidity, 3.0 < y < 3.2, our acceptance permits also the 
determination of the mean p^. of deuterons, {px)d- To obtain this results we use the deuteron 



transverse momentum spectra which were analyzed in |5^. At beam rapidities deuteron 



spectra were found to be approximately thermal with an effective (Boltzmann) temperature 
of about 80 MeV for centralities corresponding to our centrality regions 2-3, and with slightly 
lower effective temperatures of about 50-60 MeV for the most central collisions. Using the 
deuteron spectra and vi{pt) from Fig. |^ (centrality regions 1-2), one obtains values 
of {Px)d = Vi{pt)d ~ 270 MeV/c, with a systematic relative error of about 20% (coming 
from uncertainties both in the spectra and in Vi{pt)). This value of {px)d is approximately 
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twice the value of {px)p for the same centrahty and rapidity regions [0. This observation 
is in accordance with the results of Ref. []TB|. The observed absolute values of {px)/A ^ 



130 MeV/c are similar to the values {px)/A ^ 120-135 MeV/c [|T2| observed at a beam 
kinetic energy of about 1 GeV per nucleon and significantly larger than the values measured 
at still lower beam energies (cf. for instance the result of {px)/A ^ 60 MeV/c |^ measured 
at a beam energy of 0.2 GeV per nucleon). 

In Fig. ^, we compare the fiow signals of different light nuclei in the beam rapidity region 
for different centralities of the collision. The flow signal increases with mass number, flow 
of tritons and '^He is the same within error-bars. There is a rather large difference between 
proton and deuteron flow, while for nuclei with A > 2 the mass dependence is weaker and 
appears to saturate. This is very similar to the mass dependence of flow at much lower beam 
energies of 0.2 GeV per nucleon 0, where the value of (px/pt) was also analyzed in the beam 
rapidity region . Note the very high values of Vi (of the order of 0.7-0.8) for fragments with 
large transverse momenta. 

The error-bars shown in Figures |] and ^ represent statistical errors only. The systematic 
uncertainties are dominated by three sources: i) Possible particle misidentiflcation. It leads 
to relative errors in vi of approximately 5% (see above), ii) The uncertainty in the determi- 
nation of the reaction plane resolution [T^|l^, leading to a relative error in vi of the order 
of 5-10%, similar for all particle species, iii) The uncertainty in correction for flnite detector 
occupancy (see [|T^]). The accuracy of the correction itself we estimate to be of the order of 
20-30%. The occupancy correction strongly depends on the particle transverse momentum. 
The correction is maximal for the lowest pt points, where it reaches the (absolute) values 
of 0.1-0.12 (with an uncertainty of the order of 0.03-0.05). This large uncertainty in the 
occupancy corrections in the spectrometer region close to the beam limits our measurements 
at very low pt. The occupancy correction is negligible at pt >0.6-0.8 GeV/c. 



IV. DISCUSSION 



In a simple picture where the directed flow is solely due to a common collective motion 
of the matter one would expect that for light fragments {Px)a = A{px)p. At lower beam 
energies this equality was found to be approximately valid at rapidities close to the projectile 
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rapidity At lower rapidities, the ratio {px/ A) / {'Px)p was found to be increasing with 

particle mass. For A < 4 the dependence is rather significant; d{{px/A))/dy may increase 



by almost a factor of 3 from protons to He [[12| , p!3| . In our study we analyze the transverse 



momentum dependence of directed flow. Note that an equality {px)a = A{Px)p does not 
imply that Vi{pt) = Av{{pt). In the (sideward) moving thermal source model pp| , |25[ | fi(pt) 
to flrst order does not depend on the mass of the particle, and, for example, vf{pt) = Vi{pt). 
Still, in the same model {px)d = '^{Px)p, which is due to the increase of (pt) with the mass of 
the particle (assuming also a linear dependence of vi on pt)-^ Also note that in the simple 
coalescence model where in order for nucleons to coalesce one requires that they are close to 
each other only in momentum space, vf^{pt) ~ Avi{pt/A). Taking into account that "^1(^4) 
depends almost linearly on pt one arrives once more at the equality vf{pt) ~ Vi{pt). The 
experimental data (Fig. 2) show that, in the rapidity region y < 2.6, vf{pt) is indeed close 
to Vi{pt). However, at larger rapidities, we observe a significant excess of deuteron fiow in 
comparison with that of protons. It could imply that volume effects (to form a deuteron both 
nucleons should be close to each other not only in momentum space but also in configuration 
space) and/or projectile fragmentation processes become significant in deuteron production. 

In order to check if a coalescence picture including volume effects could account for our 
data we have used the coalescence model of Ref. |^ combined with the RQMD (version 



2.3) event generator [Q. This coalescence model explicitly requires the nucleons, in 
order to coalesce, to be close both in the momentum and in the configuration space. The 
closeness is defined by the cluster wave function, the parameters of which were determined 
by measured root mean square charge radii of the clusters (of the order of 1.5-2.0 fm). The 
results of our calculations of Vi{ptY'^ corresponding to the centrality region 1 are shown 
in Fig. 4. One can see that the model describes most of the features observed in the data 
(Fig. 2). The exception is, as was already observed in proton fiow |T^, that the data exhibit 



an approximately linear dependence of vi on pt at all rapidities, while the model shows a 
rather fast saturation of Vi in the projectile rapidity region. More relevant for the present 



^ More precisely {px)d/ {Px)p ~ {Td/Tp){vi/vi), where and Tp are deuteron and proton inverse 
slope constants and the ratio (vf/v^) stands for the average ratio of deuteron and proton vi{pt). 
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discussion though, the model does describe the fact that vf{pt) is very similar to fi(pt) at 
rapidities y < 2.6 and that the difference becomes large at rapidities close to beam rapidity. 
Looking into details of the nucleon freeze-out configuration space distribution (used for the 
coalescence) one can notice that, for rapidities y > 2.6, nucleons emitted in flow direction 
come from a signiflcantly narrower spatial distribution (especially along the flow direction) 
as compared to those emitted in the opposite direction. A smaller width of the spatial 
distribution means a higher probability for nucleons to coalesce. As a consequence, the 
deuteron distributions get an extra asymmetry beyond what is due to the asymmetry in 
momentum space of the nucleon distributions. This may account for the effect observed 
both experimentally and in the model. 

The difference in proton and deuteron flow as a function of rapidity could be also due to 



the presence of deuterons from projectile fragmentation. In a detailed comparison |22| of the 
coalescence model with the E877 experimental data on deuteron production it was shown 
that the model fails to reproduce the deuteron pt spectra in the projectile rapidity region, 
while describing the data fairly well at smaller rapidities. It suggests that the deuteron 
production mechanism at beam rapidity may not be (dominantly) coalescence. At present 
a quantitative description of the fragmentation process at the AGS energies is not available 
yet. However, qualitatively one could argue that deuterons from projectile fragmentation 
which are concentrated close to beam rapidity exhibit stronger flow as they are less distorted 
by thermal motion, in line with the trend observed in the data. 



V. CONCLUSION 

In summary, the directed flow of light nuclei has been measured in Au+Au collisions at 
AGS energies in the forward rapidity region. The effect has been analyzed as a function of 
particle transverse momentum for different centralities of the collision. The largest flow has 
been observed in collisions of approximately half overlapping nuclei. In such collisions and 
at rapidities close to the beam rapidity region, light nuclei exhibit very strong directed flow, 
corresponding to Vi ~ 0.7-0.8 at particle transverse momenta of about 1 GeV/c. 

The directed flow (vi) increases with the mass of the light nucleus. Coalescence model 
calculations show that the increase in vi could be accounted for by asymmetries in the 
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distributions of nucleons (forming the light nucleus) both in momentum and in configuration 
space. 
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FIGURE CAPTIONS 

1. The E877 apparatus. 

2. Transverse momentum dependence of the first moment {vi) of the deuteron (filled 
circles) and proton ||19[ (open symbols) azimuthal distributions for different particle 
rapidities and centralities of the collision. 

3. Transverse momentum dependence of directed flow fi(pt) of protons, deuterons, tri- 
tons, ^He, and "^He for different centralities of the collision. All particles are from the 
rapidity region 3.0 < y < 3.2. 

4. Transverse momentum dependence of vi{pt) of protons and deuterons in the coales- 
cence model combined with the RQMD event generator. The centrality corresponds 
to the experimental region 1. 
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FIG. 1. The E877 apparatus. 
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FIG. 2. Transverse momentum dependence of the first moment (vi) of the deuteron (filled 
circles) and proton (open symbols) azimuthal distributions for different particle rapidities and 
centralities of the collision. 
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^He, and ^He for different centralities of the collision. All particles are from the rapidity region 
3.0 < y < 3.2. 



13 



2 . 4<y<2 . 6 


2 . 8<y<3 

T 


3 . 2<y<3 . 4 ^ 

^ • deuterons 
^ o protons 


-1 

< 


> 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 



0.5 1 0.5 1 0.5 



Pt (GeV/c) 
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